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Abstract 

In this work the fohowing lepton flavor violating r and fi decays are studied: r~ fj,~ fi~ fj,'^ , 

o: 

, r — >eee+,/i ^eee+,T -^fi — >e 7 and /i ^67. We work in a supersymmetric 

(N 

scenario consisting of the minimal supersymmetric standard model particle content, extended by 
the addition of three heavy right handed Majorana neutrinos and their supersymmetric partners, 
. and where the generation of neutrino masses is done via the seesaw mechanism. Within this context, 

a signiflcant lepton flavor mixing is generated in the slepton sector due to the Yukawa neutrino 



in 
o 

o 



p. 



couplings, which is transmited from the high to the low energies via the renormalization group 
equations. This slepton mixing then generates via loops of supersymmetric particles signiflcant 



I contributions to the rates of Ij 3/j and the correlated Ij — > /j7 decays. We analize here in 

o 

full detail these rates in terms of the relevant input parameters, which are the usual minimal 



supergravity parameters and the seesaw parameters. For the Ij — > 3/j decays, a full one-loop 
QJ , analytical computation of all the contributing supersymmetric loops is presented. This completes 

and corrects previous computations in the literature. In the numerical analysis compatibility with 



^ ■ the most recent experimental upper bounds on all these r and ^ decays, with the neutrino data, and 

with the present lower bounds on the supersymmetric particle masses are required. Two typical 
scenarios with degenerate and hierarchical heavy neutrinos are considered. We will show here that 
the minimal supergravity and seesaw parameters do get important restrictions from these r and /j, 
decays in the hierarchical neutrino case. 
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I. INTRODUCTION 



The present strong evidence for lepton flavor changing neutrino oscillations in neutrino 
data implies the existence of non-zero masses for the light neutrinos, and provides the 
flrst experimental clue for physics beyond the Standard Model (SM). These oscillations also 
give an important information on the neutrino mixing angles of the Maki-Nakagawa-Sakata 
matrix (Umns) [2J- The experimentally suggested smallness of the three neutrino masses 
can be explained in a very simple and elegant way by the seesaw mechanism of neutrino 
mass generation P|. This mechanism is usually implemented by the introduction of three 
heavy right-handed (RH) Majorana neutrinos whose masses, mjvf., can be much higher than 
the SM particle masses. The smallness of the light neutrino masses, m^., appears naturally 
due to the induced large suppression by the ratio of the two very distant mass scales that 
are involved in the 3x3 seesaw mass matrices, the Majorana matrix rriM and the Dirac 
matrix m^). For instance, in the one generation case, where the seesaw model predicts 
nil, ~ iTT-'o/rnM, light neutrino masses in the 0.1-1 eV range can be generated with mo 
being of the order of the electroweak scale, v = 17A GeV, and large of the order of 
10^^ GeV. This huge separation between ruM and the electroweak scale has, however, a 
serious drawback since it leads to the well known hierarchy problem of the SM, where a tree 
level Higgs boson mass of the order of v is driven by the radiative corrections involving the 
Majorana neutrinos to very unnatural high values related to the new scale m^. 

The most elegant solution to this hierarchy problem is provided by the introduction of the 
symmetry relating fermions and bosons, called supersymmmetry (SUSY). When the seesaw 
mechanism for the neutrino mass generation is implemented in a SUSY context, the SUSY 
scalar partners of the neutrinos, i.e. the sneutrinos, do also contribute to the radiative 
corrections of the Higgs boson masses and cancel the dangerous contributions from the 
Majorana neutrinos, solving in this way the hierarchy problem of the simplest non-SUSY 
version of the seesaw mechanism. 

The best evidence of supersymmetry would be obviously the discovery of the SUSY 
particles in the present or next generation colliders. However, there are alternative ways 
to test supersymmetry which are indirect and complementary to the direct SUSY particle 
searches. These refer to the potential measurement of the SUSY particle contributions, 
via radiative corrections, to rare processes which are being explored at present and whose 
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rates are predicted to be highly suppressed in the SM. Among these processes, the Lepton 
Flavor Violating (LFV) r and fi decays are probably the most interesting ones for various 
reasons. On one hand, they get vanishing rates in the SM with massless neutrinos and 
highly suppressed rates in the SM with massive netrinos. The smallness of these rates 
in the non-SUSY version of the seesaw mechanism for neutrino mass generation is due to 
their suppression by inverse powers of the heavy scale niM- On the other hand, although 
these decays have not been seen so far in the present experiments, there are very restrictive 
upper bounds on their possible rates which imply important restrictions on the new physics 
beyond the SM. These restrictions apply even more severely to the case of softly broken 
SUSY theories with massive neutrinos and the seesaw mechanism, since these give rise to 
higher rates being suppressed by inverse powers of the SUSY breaking scale, msusY < 1 
TeV, instead of inverse powers of tjim- 

We will be devoted here in particular to the LFV r and yU decays of type Ij ^ LTand L 
3li where the present experimental upper bounds are the most restrictive ones P, M, Q, 0, Si > 
specifically, 

BR{t- e'e'e-^ 
BR{t- 

BR{t~ e~7) 
BR{fi- e~7) 

Our aim in this paper is to analize the branching ratios that can be generated for all these 
processes in the context of the SUSY-seesaw scenario with the minimal SUSY content, i.e the 
Minimal Supersymmetric Standard Model (MSSM), enhanced by the addition of three RH 
neutrinos and their corresponding SUSY partners. These LFV processes are induced by loops 
of SUSY particles which transmit the lepton flavor mixing from the slepton mass matrices 
to the observable charged lepton sector. The intergenerational mixing in the slepton sector, 
{Mj)ij, i j, is induced in turn by the radiative corrections involving the neutrino Yukawa 
couplings Yi, or, simmilarly, by the running of the soft SUSY parameters in the slepton sector, 
via the Renormalization Group Equations (RGEs), from the high energy scale, Mx > rnM, 
where the heavy Majorana neutrinos are still active, down to the electroweak scale. We will 
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assume here a Minimal Supergravity scenario (mSUGRA) with universal soft parameters 
at Mx and the breaking of the electroweak symmetry being generated radiatively. This 
scenario is also refered to in the literature as the constrained MSSM (CMSSM). 

The above LFV processes have previously been studied in the SUSY-seesaw context 
by several authors jlO|, |ll|, , under some specific assumptions for both seesaw 

parameters, m^i (or Yy, since they are related by mr) = Y,^vsmP, with tan/5 = V2/V1 being 
the ratio between the two MSSM Higgs vacuum spectation values) and niM, and for the 
mSUGRA parameters, Mq, M1/2, Aq, sign(/i) and tan/3. 

Our present study of these decay channels updates, completes and corrects the previous 
anlayses in several respects. First we include, by the first time to our knowledge, the full set 
of SUSY one-loop contributions to the Ij 3/j decays, namely, the photon, the Z boson, and 
the Higgs bosons penguin diagrams, and the box diagrams. The most complete computation 
so far of these Ij 3/j decays was done in where the contributions from the photon and 
Z boson penguin diagrams and from the box diagrams were included, but they focused on 
the particular choice of degenerate heavy Majorana neutrinos and they presented numerical 
results just for /i ^ 3e decays. We extend this previous study by including in addition the 
Higgs penguin diagrams mediated by the three neutral MSSM bosons, Hq, and Aq, and 
correct their results for the Z penguin contributions. We also extend their study in that we 
present results for the three decays, — > 3e, r — ^ 3yU and r — > 3e and consider both possible 
scenarios, degenerate and hierarchical heavy neutrinos. 

The contributions from the Higgs penguin diagrams, in the SUSY-seesaw model, were 
firstly analized in j^]. They worked in the large tan/5 limit and used the mass insertion 
approximation to account for the induced effect from the intergenerational slepton mixing 
in the SUSY contributing loops. There it was concluded that these Higgs-mediated contri- 
butions can be very relevant in the large tan f3 region, because the radiatively induced LFV 
Higgs- r-/i couplings grow as tan^/? (and, in consequence, the BR{t — > 3/i) as tan^/3), and 
also because the SUSY one-loop contributions do not decouple in these couplings. These 
large tan/5 enhacement and SUSY non-decoupling behaviour were also found in the LFV 
H,ggs boson deca... ^ Q 3 A mo.e exKaus.ive .udv of .he . ^ 3, 

and other Higgs-mediated LFV r decays, including an estimate of the Higgs contributions. 



were done in 



17| . However, these previous numerical estimates of the Higgs contributions 
to LFV r and decays were done in the context of a generic MSSM (see also ^18^), where 
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the Higgs boson mass, or equivalently m^, is an input parameter and can take small values 
of the order of 100 GeV which produces larger rates. A more recent study on the LFV 
Higgs decays has been done in 19| in the SUSY-GUT SU (5) context. We instead work here 
in the mSUGRA context where all the MSSM particle masses are quantities derived from 
the mSUGRA parameters. We will see here that this and the requirement of compatibility 
with the present experimental lower bounds for all the SUSY particle masses [2^ do indeed 
constraint the contribution from the Higgs penguins. 

In the present work we also include the predictions for the Ij li'y channels which, for the 
context we work with, are interestingly correlated with the In — » S/j rates. This correlation 
has been studied previously in the generic MSSM context in |l7| and in a similar mSUGRA 
context in Q], but in this later the dominant photon penguin approximation was used. We 
will update this comparative analysis of the Ij — > U'-f and Ij 3li rates, in the mSUGRA 
context, including the full contributions, and considering the very recent upper bounds for 
r — i> /i7 |3] and r — 67 Q]. In addition, we also require the input seesaw parameters to be 
compatible with the present neutrino data. For this comparison with the neutrino data we 
use the parametrization first introduced in for the study of the /i ^ 67 decay. 

Our final goal will be to use the SUSY contributions to all the above LFV r and fi decays 
as an efficient way to test the mSUGRA and seesaw parameters. With this goal in mind 
we will analize here the size of the branching ratios in terms of the mSUGRA and seesaw 
parameters and will explore in detail the restrictions imposed from the present experimen- 
tal bounds. We will find that for some plausible choices of the seesaw parameters, being 
compatible with neutrino data, there are indeed large excluded regions in the mSUGRA 
parameter space. 

The present work is organized as follows. In section |H] we will review the basic aspects of 
the MSSM extended with three RH neutrinos, their SUSY partners and the seesaw mecha- 
nism for neutrino mass generation. The lepton fiavor mixing in the slepton sector and in the 
mSUGRA context will be explained in section ITTTl There we also include the exact diagonal- 
ization of the sfermion mass matrices, both in the slepton and in the sneutrino sectors. The 
analytical results of the LFV Ij 3li decays will be presented in section IIVI The numerical 
results for all the LFV r and fi decays will be presented in section Finally, section IVII 
will be devoted to the conclusions. 
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II. THE MSSM EXTENDED WITH RH NEUTRINOS AND SNEUTRINOS 



In this section we briefly review the additional basic ingredients that are needed to extend 
the MSSM in order to include three right handed neutrinos, their corresponding SUSY 
partners, i.e. the sneutrinos, and the generation of neutrino masses by the seesaw mechanism. 
We follow closely the notation of refs. jl^ . Q| to describe the SUSY-seesaw scenario and 
the connection with neutrino data. For the other sectors of the MSSM we assume here the 
standard conventions as deflned, for instance, in ^,221. 



We start with the Yukawa-sector of the MSSM-seesaw that contains the three left handed 
(LH) SM neutrinos z/£ • and three extra right handed (RH) massive neutrinos i/^j, whose 
Yukawa interactions provide, after spontaneous electroweak symmetry breaking, together 
with the right handed neutrino masses, the following mass Lagrangian containing the Dirac 
and Majorana mass terms, 

-L:^ass = l(A,(4f)M'(^'''^: I +h.c., (1) 

where. 

Here m^) is the 3x3 Dirac mass matrix that is related to the 3x3 Yukawa coupling 
matrix and the MSSM Higgs vacum expectation value, < H2 >= V2 = f sin/5 with 
V = 174 GeV, by niD = < H2 >■ The other MSSM Higgs doublet gives masses to the 
charged leptons by mi = Yi < Hi >, where Yi are the Yukawa couplings of the charged 
leptons and < Hi >= vi = v cos j3. The remaining 3x3 mass matrix involved in the seesaw 
mechanism, m^, is real, non singular and symmetric, and provides the masses for the three 
RH neutrinos 

The mass matrix is a 6 x 6 complex symmetric matrix that can be diagonalized by 
a 6 X 6 unitary matrix in the following way: 

U'^Tj^ujju _ £j-v _ (liQ^g(^fi2^^^ m^^, m^^, TUn^, 771^2, 'nT'Ns) ■ (3) 

This gives 3 light Majorana neutrino mass eigenstates i>i, with masses m^^ (i=l,2,3), and 
three heavy ones Ni, with masses rriNi (i=l,2,3), which are related to the weak eigenstates 
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via, 





j^.,, , u- 1 I and I I = :r I . (4) 

The seesaw mechanism for neutrino mass generation assumes a large separation between the 
two mass scales involved in and tum matrices. More specifically, we shall assume here 
that all matrix elements of m^) are much smaller than those of itlm-, fno << f^M, and the 
predictions of the seesaw model are then given in power series of a matrix defined as, 

^ = moml}. (5) 

In particular, the previous diagonalization of the mass matrix M'^ can be solved in power 
series of S,- For simplicity, we choose to work here and in the rest of this paper, in a fiavor 
basis where the RH Majorana mass matrix, itlm-, and the charged lepton mass matrix, m^, 
are fiavor diagonal. This means that all fiavor mixing of the LH sector is included in the 
mixing matrix Umns- By working to the lowest orders of these power series expansions one 
finds, in the fiavor basis, the following neutrino 3x3 matrices. 



= -moi^ + 0{mD^^) ~ -m]:,m^}m^ (6) 
rriN = ruM + 0{mDi) — niM- 

Here, is already diagonal, but my is not yet diagonal. The rotation from this fiavor 
basis to the mass eigenstate basis is finally given by the MNS unitary matrix, Umns- Thus, 

^diag ^ ulf^sniyUMNS = diag{my^,my^,my^), (7) 

diaq i • / \ 

and the diagonalization of M'^ in eqs. (j2)) and (jHI can be performed by the following unitary 
6x6 matrix: 
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As for the Umns matrix, we use the standard parametrization given by, 



U, 



MNS 



where c,- 



v 



C12C13 



■S12C13 



■S12C23 - Ci2S23Si3e*'' C12C23 - Si2S23Sl3e"' 



S23C13 



S12S23 - Ci2C23Sl3e*'^ -C12S23 - Si2C23Sl3e*^ C23C13 



diag{l, e*", e 



(9) 



cos 6ij and Sij 



sin 6ij . 



Regarding the sneutrino sector, and because of SUSY, the introduction of three RH 
neutrinos, zz/j, leads to the addition of the three corresponding SUSY partners, ur. Thus, 
there are two complex scalar fields ul and ur per generation, as in the charged slepton 
case where there are II and Ir. The difference is that in the sneutrino sector, the seesaw 
matrix ^ is involved, as in the neutrino sector, and gives rise to a natural suppression 
of the RH sneutrino components in the relevant mass eigenstates. This fact makes the 
diagonalization procedure simpler in the sneutrino sector than in the charged slepton one. 
In order to understand properly this feature of the MSSM-seesaw model, we will illustrate in 
the following the simplest case of one generation, where this suppression already manifests. 
For this, we follow closely |23j|]. The generalization of this decoupling behaviour of the 
components to the three generations case is straightforward and we omit to show it here, 
for brevity. 

One starts by adding the new terms in the MSSM Lagrangian that involve the ur and/or 
ur. In particular, the usual MSSM soft SUSY breaking potential must be modified to include 
new mass and coupling terms for the right handed sneutrinos, which for the one generation 
case are the following. 



soft 
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^/2mw sin/3 



(10) 



where m^, and Bm are the new soft breaking parameters. These are in addition to the 
usual soft parameters of the slepton sector, mi, and Ai. The sneutrino mass terms of 
the MSSM-seesaw model can then be written in the one generation case as. 

Ml 



'~'mass 



Re{PL) Reii^R) Im{PL) Irn{pR) 



Ml 



Rc^Vr) 

Im{i'L) 
\ Im{uR) ) 
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Ml 



with, 

(m| + m|, + |m| cos 2/3 m£)(Aj, — /i cot /? ± mjv/) | ^^^^ 
moiAu - yu cot /3 ± -m^-^ + + ± 2BMmM J 
Notice that, in the sneutrino sector, there are several mass scales involved, the soft SUSY- 
breaking parameters, mi, rrij^j, Bm and A,j, the Dirac mass m/), the /x-mass parameter, the 
Z boson mass mz and the Majorana neutrino mass mM- Our basic assumption in all this 
work is that mM is much heavier than the other mass scales involved (except Mx), mu » 
mD,'mz, ^J',m,i,mj^, A^, Bm- The size of Bm has been discussed in the literature j22| and 
seems more controversial. For simplicity, we shall assume here that this is also smaller than 
mM- In this large mM limit, the diagonalization of the previous sneutrino squared mass 
matrix is simpler and leads to four mass eigenstates, two of which are light, ^2 and two 
heavy, ^2 • In the leading orders of the series expansion in powers of ^ the mass eigenstates 
and their corresponding mass eigenvalues a re g iven by (We correct in the definitions of 
and ^2 some typos with wrong signs of ref. 

el = V2{Re{uL)-^Re{uR)) ; = V2 {Im{uL) + ^Im{uR)) 

= V2{Re{uR)+^Re{i>L)) ; ^2 = (Mz^r) - e/m(z>i)) 
X 

m^i = m| + -m| cos2/? =1= 2m£)(A,^ — /iCot/3 — i?Ar)e 

m^h^ = mM ±2BMmM + m^M + 2ml, (13) 

Here we can see that the heavy states couple very weakly to the rest of particles 

of the MSSM via their 9^ component, which is highly suppresed by the small factor ^ and, 
therefore, it is a good approximation to ignore them and keep just the light states ^[2) 
which are made mainly of ul and its complex conjugate z/^. One says then that the heavy 
sneutrinos decouple from the low energy physics. 

The generalization of the previous argument to the three generations case leads to the 
conclusion that, in the seesaw limit, ^ <^ 1, the physical sneutrino eigenstates, z/^ (/5 = 1, 2, 3) 
are made mainly of the ul i states with I = e, fi, t respectively, and their corresponding 
complex conjugates. The process from the weak eigenstates to the mass eigenstates is 
simplified to the diagonalization of a 3 x 3 sneutrino mass matrix. This is to be compared 
with the more complex case of charged sleptons where the corresponding process requires 
the diagonalization of a 6 x 6 slepton mass matrix. This will be presented in the next section, 
where the most general case with lepton flavor mixing is considered. 



To end this section, we shortly comment on the parameterization that we use to make 
contact with the neutrino data. It was first introduced in [3] to study the /i — > 67 decay and 
used later by many other authors. The advantage of this parameterization is that instead of 
using as input parameters the seesaw mass matrices m^i and rriM it uses the three physical 
light neutrino masses, m^^, the three physical heavy neutrino masses, rrij^., the Umns matrix, 
and a general complex 3x3 orthogonal matrix R. With our signs and matrix conventions, 
the relation between the seesaw mass matrices and these other more physical quantities is 
given by, 

ml=im'r'''Rmt-^'l'Ut,^s (14) 

where R^R = 1 and, as we have said, m^r^ ~ "ruMi- Thus, instead of proposing directly 
possible textures for mo, or Y^, one proposes possible values for 112^-^ ,^1^^ ,rnN^ and R, and 
sets mj,^ ,m^2 y^ua and Umns to their suggested values from the experimental data. Notice 
that any hypothesis for R different from the unit matrix will lead to an additional lepton 
flavor mixing, besides the one introduced by the Umns- Notice also that the previous relation 
holds at the energy scale rriM, and to use it properly one must apply the Renormalization 
Group Equations to run the input experimental data mf"-^ and Umns from the low energies 
mw up to mjvf- Therefore, we will also include these running effects in the numerical results 
for all the branching ratios presented in this work. 

Regarding the matrix R, we will consider the following parameterization: 



^ C2C3 — C1S3 — S1S2C3 S1S3 — C1S2C3 



R 



(15) 



C2S3 C1C3 — S1S2S3 — S1C3 — C1S2S3 
y S2 S1C2 C1C2 y 

where q = cos^j, Sj = sin^j and 61, 62 and ^^3 are arbitrary complex angles. This parameter- 
ization was proposed in ref. |3l for the study of /i — 67 decays. It has also been considered 
in ref. with specific values for the 6i angles to study the implications for baryogenesis 

in the case of hierarchical neutrinos. And it has also been considered by to study the 
LFV Higgs boson decays into ijj. 

Finally, for the numerical estimates in this work, we will consider the following two 
plausible scenarios, at the low energies, being compatible with data: 
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Scenario A: with quasi-degenerate light and degenerate heavy neutrinos, 

m^, = 0.2 eV , m^, = m^, + — ^ , = m^, + — (16) 

Scenario B: with hierarchical hght and hierarchical heavy neutrinos, 



m^^ ~ eF , = y Am^^, , m^g = \/ Amlt^, (17) 
In the two above scenarios, we will fix the input low energy data to the following values. 



Amli = 0.008 eV, ^/A^ = 0.05 eV, 9u = 0,oi = 30°, 9^3 = 9atm = 45°, ^^13 = 0° 
and 6 = a = f3 = (See for instance, ref. Some results will also be presented for the 

alternative choice of small but no n- vanishing 613. 



III. GENERATION OF FLAVOR MIXING IN THE SLEPTON SECTOR 

Once the three z/^j and the three Ufi are added to the MSSM particle content, lepton fiavor 
mixing is generated in the slepton sector. This can be seen as the result of a misalignment 
between the rotations leading to the mass eigenstate basis of sleptons with respect to the 
one of leptons, which is generically present in the SUSY-seesaw models. This misalignment 
is radiatively generated from the Yukawa couplings of the Majorana neutrinos and can be 
sizable, in both, the charged slepton and sneutrino sectors. Usually, it is implemented via the 
Renormalization Group Equations (RGEs), which we take within the context of mSUGRA 
extended with three right-handed neutrinos and their SUSY partners. In consequence, we 
assume here universal soft-SUSY-breaking parameters at the large energies Mx >> ttim-, 
which must now include the corresponding parameters of the neutrino and sneutrino sectors, 
namely, 

iAi)ij = Ao(y/).„ iA,)ij = AoiY,)ij, I,] = 1, 2, 3 (18) 

Here, Mq and Aq are the usual universal soft SUSY breaking parameters in mSUGRA, 
{Yi)ij = Yi-6ij with Yi- = mi-fvi, and (Y^)ij = {rn£,)ij/v2. Notice that we have used the 3x3 
matrix form with i, j = 1, 2, 3 or equivalently i, j = e, /x, r. 
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The effects of the running from Mx down to tum on the soft mass matrices of the slepton 
sector are found then by solving the RGEs which now include the corresponding terms and 
equations for the Yukawas of the neutrinos and soft breaking parameters of the sneutrino 
sector, as they are active particles in this energy range. Below the energy scales ntMi the 
right handed neutrinos decouple and the effects of running from itlm down to the electroweak 
scale on the various parameters are obtained by solving the RGEs but now without the terms 
and equations containing the Yukawas and soft breaking neutrino parameters. The obtained 
values at the electroweak scale of the various SUSY parameters are the relevant ones in order 
to build the slepton and sneutrino mass matrices that will be presented below. 

To solve numerically the RGEs we use the Fortran code SPheno js^] that we have adapted 
to include the full flavor structure of the 3x3 soft SUSY breaking mass and trilinear cou- 
pling matrices and of the Yukawa coupling matrices. This program solves the full RGEs 
(i.e. including the commented extra equations and neutrino terms) in the two loops ap- 
proximation, computes the MSSM spectra at low energies, and uses as inputs the universal 
mSUGRA parameters, Mq, Aq, M1/2; the value of tan/5 at the electroweak scale, and the 
sign of the /i mass parameter. The value of Mx is derived from the unification condition 
for the SU{2) and U{1) couplings, gi = g2- For all the numerical analysis performed in this 
work, we have got very close values to = 2 x 10^^ GeV. The value of I fi\ is derived from 
the requirement of the proper radiative electroweak symmetry breaking. 

We present next the slepton mass matrices, relevant to low energies, that include the 
lepton mixing generated from the neutrino Yukawa couplings by the RGEs. For the charged 
slepton case and referred to the (e^, e/j, /2l, /i/?, t'l, f^j) basis, the squared mass matrix can 
be written as follows, 



/ Mlf Ml^^ Mlf Ml^^ Mll^ Wjjjl \ 
Mf,f M^f M^j,^^ Ml^^ MJ^2 

M^f M^^' M^f M^f M^l^ M^^' 
M^f Mf^f M^f M^T' Ml.l^ 

Mlf Mlf Mll'^ Mlf Mll^ Ml}^^ 

V Mlf MJtn' M^f MZ' ^..2 



(19) 
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where, 



rrif .. + 



vl {yiYi) .. + ml cos 2/? (^-^ + sin^ 9w 



m\ .. + vl ( Y^Yi ) - m| cos 2/3 sin^ dw^ij 



(20) 



The soft SUSY breaking mass matrices and trihnear couphng matrices above, m 
and A^i\ with i,j — e, fj,,T, refer to their corresponding values at the electroweak scale 
which we get with the SPheno program. 

After numerical diagonalization of the M? matrix one gets the physical slepton masses 

and the six mass eigenstates {h, , ^5)= ^ which are related to the previous weak eigenstates 

(eL,....fji)= I' by I' — R^^H, where i?^'^ is a 6 x 6 rotation matrix such that, 

^L. = R^'^MfR^^^^ = d^ag{ml..,ml). (21) 

For the sneutrino sector, the 3x3 squared mass matrix, referred to the v' — {i'cL, ^t,l) 
basis can be written as follows, 

/ m| +^m|cos2/3 —'^ —'^ ^ 



M- = 



ml 

m~ m~ +^m7 cos 2/3 



L,eT 

L,fJ,T 



\ 



m 



L,Te 



m 



L.Tfl 



(22) 



where m| ^. are the same as in the previous charged slepton squared mass matrix. After 
diagonalization of the M| matrix one gets the relevant physical sneutrino masses and eigen- 
states, i/p {(3 — 1,2,3) which are related to the previous states D'^ by the corresponding 3x3 
rotation matrix, v' — R^'^^u, and is such that, 

Ml^^^ = R^^^MjR^^^'^ = dtag{ml,ml,ml). (23) 

Finally, in order to illustrate later the size of the misalignment effects in the slepton sector 
we define the following dimesionless parameters, 

(24) 
(25) 
(26) 















Xij _ 










w? 






"RR — 
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where 




) 



1/6 



(27) 



is an average slepton squared mass. These parameters have also been considered by other 
authors in a more model independent approach and with the purpose of getting bounds from 
experimental data. Some of these bounds can be found in [3, 1^ • 

For all the numerical results presented in this paper, we will set values for all the following 
input parameters and physical quantities: 

• mSUGRA parameters: Mq, M1/2, Aq, sign(/i) and tan/3. 

• seesaw parameters: rriNi, 1TIN2, n^'N-j, and R (or equivantly 9i, 6*2, 6*3). 

• physical quantities: m,^^, m^^, rriu^, Umns 

IV. ANALYTICAL RESULTS FOR THE r ^ I'ljlf DECAYS 

J ILL 

In this section we present the analytical results for the LFV r and yU decays into three 
leptons with the same flavor, within the mSUGRA-seesaw context that we have presented 
in the previous sections. We perform a complete one-loop computation of the r and [i decay 
widths for all the three possible channels, — ^ r — > e^e^e+ and /i — > e^e~e"'", and 

include all the contributing SUSY loops. We present each contribution separately, 7-penguin, 
Z-penguin, Higgs-penguin and boxes. The contributions from the Higgs-penguin diagrams 
are, to our knowledge, computed exactly by the first time here. We have also reviewed the 
analytical results in [lO| and correct their results for the Z-penguin contributions. Notice 
that we make the computation in the physical mass eigenstate basis. That is, we consider the 
one-loop contributions from charged sleptons, Ix (-^ = 1, ••,6), sneutrinos Vx {X = 1,2,3), 
charginos Xa ~ 1)2), and neutralinos Xa {A = 1,..,4). In all this section we follow 
closely the notation and way of presentation of The interactions in the physical mass 
eigenstate basis that are needed for this computation are collected in the form of Feynman 
rules in Appendix 1X1 

First, we define the amplitudes for the lj{p) {Pi)^i^ {P2)^t {Ps) decays as the sum of 
the various contributions, 



T(/7 ^ /r/r/+) = T 



7— penguin 



+ Tz. 



—penguin 



+ Th 



—penguin 



~l~ -^boxes' 



(28) 
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In the following we present the results for these contributions in terms of some convenient 
form factors. 



A. The 7-penguin contributions 

The diagrams where a photon is exchanged are called 7-penguin diagrams and are shown 
in fig. lU The result for the 7-penguin amplitude contributing to the Ij KK^t decays is 
usually written as, 

T^-penguin = Hvi) [q'lMiPL + A^Pr) + im/^.a^,g" [A^Pl + A^Pr)] Uj{p) 

_ 

X -^Ui{p2h^Vi{p3) - {pi ^ P2) (29) 

where q is the photon momentum and e is the electric charge. The photon-penguin ampli- 
tude has two contributions in the MSSM-seesaw from the chargino and neutralino sectors 
respectively, as can be seen in the structure of the form factors, 

j^L,R ^ j^(n)L.R ^ J^{c)L,R^ a = 1, 2 (30) 

The neutralino contributions are given by 

{n)L _ 1 ^jR ^rR* 1 2 - 9xax + I8x\j^ ~ llx\ + 6x^^ log xax 



' 32772 ^2 



.rL .jL* 1 - QXAX + ^Xax + 2xjtx - 6X^X log 3^ AX 

<6{1-xax) 



(31) 



, .tR mR* ^ ~ ^^^^ + ^^^^ + ^^^^ ~ ^^^^ log 

+ ^^lAX^^iAX' 



(32) 

^in)R ^ (33) 



, ArL nrR* ^ - ^ AX + AX log X AX 

+ '■^iAx^^jAx— ::; ^ 

m, (1 - Xax) 



where xax = ''^-o /'^r • On the other hand, the chargino contributions are 

Xa' Ix " 

r-iR r^R* 1 16 — 45xax + 36x^^ — 7x^4 + 6(2 — 3xax) logx^x 
' " 576.2 (1_^^^)4 

(34) 
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(G5) 



kiP'i) hip) *{ 




FIG. 1: 7-penguin diagrams contributing to the 1^ Z- 1^ if decay 



327r2 ml 



6(1 -xax)^ 



rR* mu2 + ^XAx-Qx\x + x\x + QxAx\ogXAx 

+ ^iAX^jAX' 



+ ^iAX^jAX 
a(c)R _ a{c)L\ 



R* 



6 (1 - xax) 

-3 + Ax AX - Xax - 2 iogXAX 
(1 - XAxf 



(35) 
(36) 



where xax = ml-fml . Notice that in both neutrahno and chargino contributions a sum- 



Xa' 



mation over the indices A and X is understood. Notice also that we have not neglected any 
of the fermion masses. If we neglect these masses in the previous formulas we get the same 
result as in jl^. The expressions for the and C couplings are given in the Appendix 1X1 
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B. The Z-penguin contributions 



The diagrams where a Z boson is exchanged are called the Z-penguin diagrams and are 
shown in fig. |21 The amplitude in this case is 

^z-penguin = i^i(Pi) (FlPl + FrPr)] Uj{p) 



where, as before, -Fl(_r) = ^ij^r) + ^l{r)- '^^^ 
following: 



'ife) - {Pi ^ P2) (37) 
expressions for these form factors are the 



1 



+ N^AxNfXy 



2(5xy<^24(m|o^,m?^,m?^) 



+ N^AxNfXx 



^ L 



(38) 
(39) 



167r2 



2(5xy ^"24 {mi- , m?^ , m?^ ) 



"T ^iAX^jAX 



_ p(c) 
^R — 



L^R 



(40) 
(41) 



Notice that all the loop functions are evaluated at zero external momenta which i 
good approximation in these decays. That is. 



C{ml,ml,ml) = C(0, 0, m?, m^, m^) 



is a very 

(42) 
(43) 



The expressions for the couplings are collected in Appendix^ and the loop functions 31 1 
are given in the Appendix El Notice that our result for the Z-penguin contributions differs 
significantly from the result in ^L^. In fact, these authors did not consider all the diagrams 
in these Z-penguin contributions, which we think is not justified. 



C. The box contributions 



The box-type diagrams are shown in fig. El We have computed these diagrams and found 
a result in agreement with [10]. The amplitude for these box contributions can be written 
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Tboxes = e''B^[Ui{pi){YPL)Uj{p)][Ui{p2){'y^,PL)Vi{ps)] 

+ e^B^ [uM (rPR) Uj{p)] [uM {i^Pr) vM\ 

+ e^B^ (T'^Pi) u,{p)] [u,{p2) [i^^Pr] t^.fe)] - {Pi ^ P2)] 

+ e'B^ {[«,(pi) {^Pr) u,{p)] [Ui{p2) {h^Pl) Viip^)] - (pi ^ P2)} 

+ e'^B^ {[Ui{Pl)PLUj{p)\ [Ui{p2)PLVi{pz)\ - (Pl ^P2)} 

+ e^Bf {[Ui{pi)PRUj{p)] [Ui{p2)PRVi{p^)] - {pi ^P2)} 

+ e^Sf {[Mi(pi) {a^^PL) Ujip)] [u^{p2) (a'^'-PL) v^{p3)] - (pi ^ P2)} 
+ e^Sf ((7^,P^ii,(p))] [iZ,(p2) ((j'^'^Pi?) Viips)] - (pi ^ P2)} 
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FIG. 3: Box-type diagrams contributing to the I- —>■ L L decay 



where 

The different neutrahno contributions are, 



1,...,4 



g2^(n)L 



2 r>{n)L 



e'B. 



2 r>in)L 



where 



1 



167r2 



1 



167r2 



^ -HAY^^ jAX'-HBX'-HBY 2 



Xa Xb iAY^ \jAX^HBX^HBY 



^ ^^iBY^^iBX^^jAX^^iAY + ^ iBY^^ iBX^^ jAX^^ iAY 



167r2 

1 



Dom^om^.^Nr^yNfXxNrBxNfBY + ^^x^^x?,^"ByA^"BxA^"AA^My 



(45) 



(46) 



(47) 



R* atR* 



^m^om^o^NfXxN^'^YN^BYNl'BX 



1,...,4 



(48) 
(49) 
(50) 



Do 
Do 



7^/ 2 2 2 2 ^ 



(51) 
(52) 
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The chargino contributions read, 



167r2 

1 

167r2 







^0 f^R ^R* ^R ^R* 
~^'~'iAY^jAX^iBX^iBY 



(53) 



'^^lAY^jAX^iBX^iBY ^ Xa X~b iAY^ jAX^iBX^iBY 



fL* 



R* r-iR /^L* 



Dom^- m - - Cfj^Y^flx^iBxCiBY 



a \l^R 



1,...,4 



(54) 
(55) 
(56) 
(57) 



where 



^ ^ Xa' Xfl' 



Xa Xb 



(58) 
(59) 



D. The Higgs-penguin contributions 



The diagrams where a Higgs boson is exchanged are called the Higgs-penguin diagrams. 
These are shown in fig. 0] and have been computed here by the first time. These are usually 
not considered in the literature. In particular, in the most complete study so far of 10| 
these Higgs-penguin diagrams were not included. However, they are expected to be relevant 
at large tan/3 We will therefore include them here. Specifically, we include the contri- 
butions from the three neutral MSSM Higgs bosons, /iq, Hq and Aq and consider all SUSY 
loops. 

In this case, the amplitude can be written as. 



Higgs 



2 tdL 



2,Higgs 



+ e'B. 



2 uR 



2,Higgs 



{[Ui{pi) {YPl) Uj{p)] [u,{p2) (i^Pr) V.,{P3)] - {pi ^ P2)} 
{[ui{pi) {i^Pr) Uj{p)] [ui{p2) {i^Pl) Vi{pi)] - {pi ^ P2)} 



+ e 53,Higgs{K(Pl)^-L«j(p)] [Ui{.P2)PLVi{p^)] - (pi 



+ e^^Jffiggs {[UiiPl)PRUj{p)] [Ui{p2)PRVi{p3)] - (pi 



P2)} 
'P2)} 



(60) 



where 



T^L,R _ r^{n)L,R rj{c)L,R _ 9 o 

-°a,Higgs ~ -°a,Higgs -°a,Higgs " — ^, -J 



(61) 
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h{Pi) 
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kiPs) 



(H3) 



(H4) 




FIG. 4: Higgs-penguin diagrams contributing to the Ij — > Zj if decay. Here Hp{p = 1, 2, 3) 



The first term represents the neutrahno contribution, which we find to be 

3 



6 -°2,Higgs 



2D(n)i 
e -°3,Higgs 



= E 



B 



(n)R 



p=l 
3 



'2) ml ^'"^'^ 



^2 ^L,n'^L,i 



B 



(n)L 



a,Higgs a,Higgs 



a = 2,3 



(62) 

(63) 
(64) 
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where Hp{p = 1, 2, 3) = H^, and 



H 



D/2 2\| 2^/2 2 2\| 2r^/2 2 2 \ 



^^iAX^R,AB^^ jBX 



+ (Cll - C'12) ("^[^ , "^|o^ , "^|o^ ) 



R* 
jBX 



(P) t\tL* 
R,AB^^jBX 



+ "^xi"^X°s'^o("^L ' ^hB^^^^AxDi^AB^fSx 

(p)'' r ™ /'/^ ^/'™2 _2 _2 ^ ArH* 



-AT, 



+ 



nip) 



2 2 

mf — mf 



jAX 



+ 



:(p) 



5 



2 2 

mf — mf 



-m^ 5i(m|o^, m?^)iV^^iVf;^ + mz^.m^o^5o(m|o , m?^)iV^^iV/;^ 



(p) 



m,,m,^,5i(m|c, , m?_^)iV^^^iV/;^ + mz,m^^So(m|^, m?^)iV^^^iVf;;,] } (65) 
p= 1,2,3 (66) 



if 



(p) 



The values of the couphngs are given again in Appendix El and the loop functions in Ap- 
pendix ^ Correspondingly, the result for the chargino contribution is given by, 

3 



2 ^(c)L 
^ ^2,Higgs 



2o(c)i 
^ -°3,Higgs 



p=l 
3 



/^(P) qiP) 
ml ^'-^^'^ 



p=i 



B 



{c)R 
a,Higgs 



B 



{c)L 
a,Higgs 



L'^R 



a = 2,3 



(67) 

(68) 
(69) 
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where H^ur) c ^^"^ obtained from the previous -^^1%) „ by replacing everywhere 



[ - 








x' - 








Dl(r) - 





Again the values of the couplings and the loop functions are given in Appendices El and iBl 
respectively. 



E. L L L if decay width 



The decay width for I- I- I- If can be written in terms of the form factors given in the 
previous sections as [l^ : 



5127r3 



|Af I V |Af 1^ - 2 {A^A^* + A^A^* + h.c.) 



i|2 



16 , rriL 
3 mi^ 



22\ 



3 / 



+ 



+ 



B^\' + \B^\' 



+ 



B. 



R 



1 

24 
1 

2 
1 



4" 



Bi(B^* + Bi^BT + h.c 



Bi 



R r>R* 



6 fii 



\B 



R\ 



+ - yA{B^* + 5f * + A{B^* + 5f * + /i.e. 
- ^ [a^B^* + A^5f * + A^^^* + A^]3^* + /i.e. 
+ ^ {2 {\Fll\' + + IF^hI' + \Frl? 



+ 2 (Af + F^^ + /..c.) + {A[Fl^ + F^^ + /..c.) 
- 4 (A«F£^ + A^F^^ + h.c) - 2 (A^F^^ + A^Fl^ + /..c.) }] 



(70) 
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where 



Fll — 


^L^L 




Frr = 


FllIl^r 


Flr = 




g"^ sin^ Oy/Tn^ 


Frl = 


FlrIl^r 



(71) 
(72) 

(73) 
(74) 

Notice that we have put the Higgs contributions together with the box ones in order to 
follow closely the way of presentation of 

B^'^ = B^'^ + B^;^,^^, (75) 
Notice that we have corrected the result in ref . for the term that goes with ( | A2 | + | A2 | ) . 



V. NUMERICAL RESULTS FOR THE LEV BRANCHING RATIOS 

We present in this section the numerical results for all the branching ratios of LFV r and 
/i decays in the context of the mSUGRA-seesaw scenario that has been introduced in the 
previous sections. We focus on the following LFV decays, t" e~e~e+ and 

yU~ —>■ e~e~e~^, and the radiative decays A* 7) ''"^ and /i~ 6^7. The reason to 

consider these radiative decays together with the decays into three leptons is that there are 
insteresting correlations among them that provide additional information in testing SUSY. 
Specifically, we show in this section the correlations between the ratios of fi^fi^fi~^ 
and T~ — s> yU~7; between r~ —* e~e~e^ and r~ — > e~7; and between /i~ e~e~e~^ and 
fi~ —>■ e~7. For the numerical estimates of the radiative decays we use the formula of jlol |. 
which is given in terms of the A2 ' as, 

nij - lii) = i^K(i^2 r + (77) 

•' iovr ^ 

but we use our expressions for the form factors in eqs. lp?^ . (jHS)), and that include the 
lepton mass contributions. We explore here in full detail the size of the SUSY contributions 
to all these LFV Ij 3/j and Ij — > /j7 decays as a function of all the mSUGRA parameters. 
Mo, Ml/2, ^0, tan/5 and sign(/i) and the seesaw parameters ttiat^, i = 1,2,3 and R or. 



24 



equivalently, 6i, 62 and 63. In all this numerical analysis we require compatibility with 
the neutrino data and with the present upper experimental bounds for all these branching 
ratios 0,0,0,18,01, as given explicitely in the introduction. We also demand the complete 
set of SUSY particle masses, which we derive with the SPheno program, to be above the 
present experimental lower bounds [2^. The numerical values of the total r and /i widths 
(lifetimes) are taken from ^J]- We show first the results for the scenario A with quasi- 
degenerate light and degenerate heavy neutrinos and next the most interesting scenario B 
with hierarchical light and hierarchical heavy neutrinos. 



A. Degenerate case 

We show in figs. El through |H1 the numerical results of the branching ratios for the LFV 
r and /i decays in scenario A with degenerate heavy neutrinos of mass m^r. We show our 
predictions for the three channels, r~ ^^fi^fi^, r~ e~e~e^ and /i~ — e~e~e^, and 
similarly, for the comparison with the leptonic radiative decays, Ij k'-f, we also show in 
the plots the correlated decay, — > fi^j, e~7 and e~7, respectively. 

The results of the branching ratios for the fi^fi^fi^ and fi~'~f decays as 

a function of tan/3 are illustrated in fig. In these plots we set m^v = 10^^ GeV and 
assume the matrix R to be real. Notice that in the degenerate case with real R these 
LFV ratios do not depend on the particular choice for R. This can be easily understood 
because the dependence on R drops in the relevant factor, (Y*Yj)ij, appearing in the 
dominant slepton mixing, and due to the property R^R = 1. From this figure we also 
see that the predicted rates for both channels are well below their respective experimental 
upper bounds for all tan/5 values, eventhough the total rates grow fast with tan/3. We also 
see clearly the mentioned correlation between the r~ ^~ and r~ /i~7 rates. 
In fact, this correlation is an inmediate consequence of the dominance of the 7-penguin 
contributions which clearly governs the size of the r~ fi~fi~fi~^ rates. This dominance 
is illustrated in fig. El (a), where the various contributions are shown separately. In fact, 
the contributions from the 7-penguin diagrams are almost undistinguishable from the total 
rates for all tan/5 values. For low tan/5 values the next dominant contribution is from the 
Z-penguin diagrams, but this is still more than one order of magnitude smaller than the 
7-penguin contribution. The contributions from the box diagrams are even smaller. We 
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FIG. 5: Dependence of LVF r decays with tan /3 in scenario A with degenerate heavy neutrinos and 
real R, for mjv = IQ-*^^ GeV. (a) Upper panel, BR{t — jx" ix~ ix^) and its different contributions, (b) 
lower-left panel, BR{t — > /X7) and (c) lower-right panel, l^'j^L lr rr\- ^he other input parameters 
are, Mq = 400 GeV, M1/2 = 300 GeV, Aq = and sign(//) > 0. 

also learn that the Z and boxes contributions do not depend significantly on tan/3, while 
the photon contribution goes approximately as (tan/3)^ at large tan/3. In this large tan/3 
region it is interesting to note that the total Higgs contribution becomes larger than the 
Z contribution and the boxes, due to the fact that it grows approximately as (tan /3)^. In 
this total Higgs contribution the dominant penguins are those with Hq and Aq exchanged, 
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which are several orders of magnitude larger than the /lo-penguin contribution. However, in 
spite of this huge enhacement of the total Higgs contribution occurring at large tan jS, its 
relative size as compared to the photon-penguin contribution is still negligible. For instance, 
for the values set in this figure of Mq = 400 GeV, M1/2 = 300 GeV, Aq = 0, sign/x > 
and m^v = 10^*^ GeV, the Higgs contribution is still four orders of magnitude smaller than 
the photon-penguin contribution at tan/5 = 50. This set of values give rise to the following 
MSSM spectrum (we just specify here the relevant sectors). 



mr 

n 


= 247 GeV 


m^o 

Al 


= 121 GeV 




= 114 GeV 


mr 


= 397 GeV 


m^o 


= 232 GeV 


irtijo 


= 457 GeV 




= 413 GeV 


m^o 


= 484 GeV 




= 457 GeV 




= 416 GeV 




= 493 GeV 


niy^ 


= 351 GeV 


mr 


= 417 GeV 




= 232 GeV 




= 409 GeV 


mr 

'6 


= 419 GeV 




= 495 GeV 




= 410 GeV. 



We have checked that other choices of parameters, specially lower Mq and M1/2 lead to 
larger contributions from the Higgs penguins, since one gets ligther SUSY spectra and more 
importantly lighter Hq and Aq bosons. However, the present experimental lower bounds on 
the MSSM particle masses, do not allow to decrease much these Mq and M1/2 values, so 
that in this mSUGRA context the relevant rriHo, ^"^^ masses can never get low enough 
values such that their corresponding Higgs-penguin contributions be competitive with the 7- 
penguin ones. From this figure we conclude then that the leading 7-penguin approximation 
works extremely well, for all tan (3 values. In this approximation one gets. 



BR{lj 

BR{lj ^ la) 




(78) 



which leads to the approximate values of ^ and for (Ijli) = (r/i), (re) and (/ie), 
respectively. As will be seen later it also works extremely well in the other channels. These 
nearly constant values of the ratios of branching ratios will be showing along this work. 
Obviously, if these ratios could be measured they could provide interesting information. 

In fig. ^c) we have included our predictions for \Sfl\, {Sf'j^l and, as defined in 

eqs. ()24j) . (f^H|) and respectively, as a function of tan/3. These are the flavor changing 
parameters that are the relevant ones for the r decays having /i in the final state. It is also 
interesting to compare them with the predictions in the leading logarithmic approximation 
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where the generated mixing in the off-diagonal terms {i ^ jihj = 1,2,3), through the 
running from Mx down to mju, is given by 

(Am|)., = -1-^{?,mI + AI){Y:LYJ\, 

{Ami),, = 0;Lki = log ( ^) 5^. (79) 

and, in consequence, it predicts the hierachy, \5'^ii\ > > 

As expected from the leading- log approximation, we see in fig. Etc) that \Sj^i^\ is much 
larger than {Sf'^l and However, we get larger than and it can be indeed 

two orders of magnitude larger than {Sf'j^l at large tan/5. It is clear that, at least for our 
choice here of = 0, the leading-log approximation does not fully work. We also learn 
from this figure that the size of the mixing is always small in the degenerate case, being the 
largest about 3 x 10^^. 

We next comment on the relevance of the choice for the mjv values. In fig. IHlwe have 
illustrated the r and r ^ yU7 branching ratios as a function of m^v for degenerate 

heavy neutrinos and tan/? = 50. The explored range in m^v is from 10^ GeV up to 10^^ 
GeV which is favorable for baryogenesis. Both rates have the same behaviour with rriN 
which corresponds approximately to BR{t /i^yU^/i+), BR{t — > /i7) oc Im^r log(m7v)p. 
As before, these two predicted branching ratios are well bellow their experimental upper 
bounds, even at the largest m^r value of 10^^ GeV. In the last plot of fig. IHlwe inlude the 
dependence of \^'rr\ on tjin which clearly show a correlated behaviour with the 

previous plots. Again, \5f^i\ is the dominant one reaching values up to about 3 x 10~^, and 
(5^^ is larger than 

For completeness, we also include the results of the other four LFV r and /i decays in fig. [3 
where the predictions are shown as a function of tan/5. These behaviours are very similar to 
those in BR{t^ —>■ ^^fi^fi^) and BR{r fi'j) decays correspondingly. The main difference 
is in the lower plots, where now l^]^^^^"*] is larger than | ^/^j^;^'^'' | • The maximum reached values 
are very small in this case, ~ 5 x 10"^. We see again that the leading 7-penguin 

approximation works extremely well for these channels, and the previously mentioned values 
of the ratios of branching ratios give a pretty good answer. We also find that the rates for all 
these four decays are well below their corresponding experimental bounds, in the degenerate 
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FIG. 6: Dependence of LFV r decays with ttt-at in scenario A with degenerate heavy neutrinos and 
real R, for tan/3 = 50. (a) Upper panel, BR{t — > /x^/i^/x^), (b) lower-left panel, BR{t fj,j) and 
(c) lower-right panel, I^'ll LR rr\- '^^^ other input parameters are, Mq = 400 GeV, M1/2 = 300 
GeV, ^0 = and sign(^) > 0. 



case, for all the explored values of taii/S and rriN. 

To end up the study of the degenerate case, we have also explored the dependence of 
the largest ratios BR{t~ fi~^~ii'^) and BR{t~ /i~7) with the mSUGRA parameters 
Mo and Mi/2- These results are shown in figlHl We see clearly a similar behaviour in the 
two channels and their rates decrease as expected when increasing the soft SUSY breaking 
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FIG. 7: Dependence of LFV r and jj. decays with tan/3 in scenario A with degenerate heavy 
neutrinos and real R, for mjv = 10^^ GeV. (a) Upper-left panel, BR{t e~e'e^), (b) upper-right 
panel, BR^jj, — e~e~e+), (c) middle-left panel, BR{t — 67), (d) middle-right panel, BR^jjl — 67), 
(e) lower-left panel, lr rr\ and (f) lower-right panel, The other input parameters 

are, Mq = 400 GeV, M1/2 = 300 GeV, = and sign(/x) > 0. 

mass parameters. This implies that for large enough values of Mq or M1/2 the branching 
ratios are considerably suppresed, due to the decoupling of the heavy SUSY particles in the 
dominant loops which are common to both observables. Thus, looking at these plots we can 
obviously conclude that the hghter the SUSY spectrum is, the larger branching ratios we 
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FIG. 8: Dependence of BR{t~ — ji' ji~ jJL^) and BR{t — /U7) with Mq and M1/2 in scenario A 
with degenerate heavy neutrinos and real R, for m^r = 10^^ GeV and tan/3 = 50. (a) Upper- 
left panel, BR{t fjrfj~fi'^) as a function of Mq for M1/2 = 100 GeV, (b) upper-right panel, 
BR{t —>■ nj) as a function of Mo for M1/2 = 100 GeV, (c) lower-left panel, BR{t — ix~ ^~ ^'^) as 
a function of M1/2 for Mo = 100 GeV, (d) lower-right panel, BR{t — /X7) as a function of M1/2 
for Mo = 100 GeV. In all the plots we take Aq = Q and sign(/Lt) > 0. 

get. However, as already said, the more interesting region of low Mo and/or M1/2 values, 
being close to 100 GeV, is not allowed by the present experimental lower bounds on the 
MSSM particle masses. 

In summary, in the case of degenerate heavy neutrinos, we get LFV r and ^ decay rates 
which are still below their present experimental upper bounds, for all the explored values 
of the seesaw and mSUGRA parameters, which have been required to provide a full MSSM 
spectrum with masses being compatible with the present experimental bounds. 
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B. Hierarchical case 



We next present the results for hierarchical neutrinos, scenario B, which are much more 
promissing. In this case the choice for R is very relevant. The results for the general complex 
R case and for the particular mass hierarchy {miy^,miy^,miy^) = (10^,2 x 10^,10^^) GeV, 
are shown in figs. lUl through IT^ This particular choice for the heavy neutrino masses seems 
to generate a proper rate for baryogenesis via leptogenesis in the hierarchical case Q]. We 
will later explore other choices as well. 

From these figures we first confirm that the LFV r and fi decay rates are much larger 
in the hierarchical case than in the degenerate one. This is true even for the case of real 
R, which corresponds in our plots to the predictions at arg(6'i) = arg(6'2) = arg(6'3) = 0. 
Furthermore, we get severe restrictions on the maximum allowed decay rates coming from 
the experimental upper bounds. 

The predictions for BR{t~ — > /i~yU~/i+) and BR{t fi'j) as a function of 16*21 are 
depicted in fig El Here 9i and 6*3 are set to zero, and arg(6^2) = 7r/4. From now on the 
arguments of 61, 62 and 6*3 are written in radians. The other parameters are set to tan (3 = 50, 
Mo = 400 GeV, M1/2 = 300 GeV, Aq = and sign(^) > 0. In fig. Efa) we show separately 
the various contributions to BR{t^ — ^ The dominant one is again the photon- 

penguin contribution (which is undistinguisible from the total in this figure) and the others 
are several orders of magnitude smaller. We also see that the relative size of the subdominant 
contributions have changed respect to the previously studied degenerate case. Now the Higgs 
contribution is larger than the boxes one and this is larger than the Z one. This is so because 
the largest tan/3 = 50 value has been set. All the rates for fi^fi^fi^ in this plot are 

within the allowed range by the experimental bound, which is placed just at the upper line 
of the rectangle. In contrast, one can see in fig Efb) that, for the chosen mSUGRA and 
seesaw parameters, the predicted BR{t —>■ fi'-f) are clearly above the experimental bound. 
The dependence of \Sj^LLRRR\ with \92\ is shown in fig. Efc). We see that l^^^.! can reach 
very large values, up to 0.4, for \92\ = 3 and arg(6'2) = vr/4. We have checked that this 
particular choice of 62 = 3e*'^/^ gives rise to large neutrino Yukawa matrix elements \Y^^\ 
and \Yj^^\ of the order of 1, which are the responsible for this large mixing in the slepton 
sector. 

It is also interesting to compare the MSSM spectrum for this hierarchical case with the 
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FIG. 9: (a) Upper panel: Dependence of BR{t~ fi~fi~fi'^) with 16*21, (b) lower-left panel: 
Dependence of BR{t — > /X7) with \62\, (c) lower right panel: \^'l'llrrr\ with 16*21. The hori- 
zontal line is the upper experimental bound. All panels are in scenario B, for arg(02) = 7r/4, 
{mN^,mN^,mN,,) = (10^2 x 10^ lO^^) GeV, 6*1 = 6*3 = 0, tan/3 = 50, Mq = 400 GeV, M1/2 = 300 
GeV, ^0 = and sign(/i) > 0. 

previous degenerate case. For the input values of fig. IHlbut with 6*2 set to the extreme value 
62 = 2.8e*^ we get the following masses, 
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It is obvious that the complex R affects significantly the predictions of the MSSM masses, 
specially in the slepton sector. In general, the slepton mixing generated by the complex 
6'j, lower the lightest charged slepton and the lightest sneutrino masses and increases the 
heaviest charged slepton and sneutrino masses. 

In fig. nil we show the predictions of BR{lJ l^l^l^) and BR{lj lij) as functions 
of \62\, for all the channels and for the different values of arg(6'2) = 0, vr/lO, 7r/8, 7r/6, 7r/4. 
In all these plots we set again tan/5 = 50, Mq = 400 GeV, M1/2 = 300 GeV, Aq = 0, 
sign(/i) > and {mjq^,mjq^,mjq.^) = (10^,2 x 10^,10^^) GeV. The upper lines correspond 
to arg(^^2) = it/4: and the lower ones to arg(6'2) = 0. These lower lines are therefore the 
corresponding predictions for real R. It is clear that all the branching ratios have a soft 
behaviour with \62\ except for the case of real 62 where appears a narrow dip in each plot. 
In this fig. El we see that all the rates obtained are below their experimental upper bounds, 
except for the processes r ^ /i7 and /i 67, where the predicted rates for complex 62 with 
large \62\ are clearly above the allowed region. The most restrictive channel in this case is 
T ^ ^'-f where compatibility with data occurs just for real 62 and for complex 62 but with 
1^2! values near the region of the narrow dip. We also see that the rates for BR{fi 3e) 
enter in conflict with experiment at the upper corner of large \62\ and large arg(^2) = tt/4:. 

Even more interesting are the predictions for BR{lJ KK^t) ^^'^ BR{lj li'y) as 
functions of |^i|, due to the large values of the relevant entries of the Y^, coupling matrix, 
which are illustrated in fig. CH Concretely, \Yj}^\ can be as large as ~ 0.2 for 16*11 ~ 2.5 and 
arg (^1) = 7r/4, and \Y^^\ and \Y/^^\ are in the range 0.1 — 1 for all studied complex 61 values. 
The results for BR{lJ l^l^lf) and BR{lj —>■ li'j) as functions of |^i|, for different values of 
arg (^1), are illustrated in fig. El Here 62 and 6^ are set to zero. The same set of mSUGRA 
parameters and heavy neutrino masses as in fig. El are taken for comparison. We see clearly 
that the restrictions are more severe in this case than in the previous one. In fact, all the 
rates cross the horizontal lines of the experimental bounds except for BR{r^ fi^fj,~fi'^) 
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FIG. 10: Dependence of LVF r and /U decays with 1^2 1 in scenario B with hierarchical heavy 
neutrinos and complex R, for arg(^2) = 0, tt/IO, 7r/8, 7r/6, 7r/4 in radians (lower to upper lines), 
{mNi,mN2,mN3) = (10^,2 x 10^, 10^'') GeV, 61 = 63 = 0, tan/3 = 50, Mq = 400 GeV, M1/2 = 300 
GeV, ^0 = and sign(/x) > 0. The horizontal lines are the upper experimental bounds. 

and BR{t~ e~e~e'^). The most restrictive channel is now the ^ decay. More 
specifically, we sec that all the points in the plot of BR{fi 67), except for the particular 
values = and real 61 at the dip, are excluded by the experimental upper bound. Also 
the predictions for BR{ii 3e) are mostly excluded, except again for the region close to 
zero and the dip. Notice that the qualitative behaviour of these all branching ratios with |^i| 
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FIG. 12: Dependence of LFV r and /i decays with \9i\ in scenario B with hierarchical heavy 
neutrinos and complex R, for arg(^i) = 0, tt/IO, 7r/8, 7r/6, 7r/4 in radians (lower to upper lines), 
(mjVi,mAr2,mjv3) = (10^2 x 10^ 10^'') GeV, 62 = 03 = 0, tan/3 = 50, Mq = 400 GeV, M1/2 = 300 
GeV and Aq = 0. The horizontal lines are the upper experimental bounds. 



a lighter MSSM spectrum and and in consequence to higher rates. For comparison with the 
previous cases, we include below the predicted masses of the relevant MSSM particles, for 
the particular value 6*1 = 2.8e*5, 
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Notice that the hghtest slepton, neutrahno, chargino and Higgs boson have masses close to 
their experimental lower bounds. 

We conclude from this fig. Uni that the predictions for BR{fi 67) and BR{jj, 3e) 
are totally excluded by present data and the predictions for BR{t —>■ fi'-f) are practically 
excluded, with the exception of the two narrow dips. The predictions for BR{t — > 67) 
get severe restrictions for complex 9i with large |6'i| and/or large arg(^i), and the rates for 
BR{t 3/i) start being sensitive to the present experimental bounds for large complex 9i 
values in the upper corner of the plot. 

We have also explored the dependence with the complex ^3 angle, and it turns out that 
the predictions for all rates are nearly constant with this angle. For instance, for tan (3 = 50, 
Mo = 400 GeV, M1/2 = 300 GeV, Ao = and sign(/i) > 0, we get BR{t = 2.6 x 10^^°, 

BR{t 3e) = 8.8 x 10"^^ BR{fx 3e) = 1.8 x 10~l^ BR{t 7x7) = 9.1 x 10"^ 
BR{t — i> 67) = 7.8 X 10"^^ and BR{jj, 67) = 2.6 x 10"^^. In this case only the prediction 
for BR{t fij) is in confiict with the experiment. 

The dependence of BR{t^ fi^fi^fi^) and BR{t with the mSUGRA parameters 

Mq and Ml/2 is illustrated in fig. We see a similar behaviour as in the degenerate case, 
where a suppresion of the branching ratios occurs for large values of Mq and/or M1/2. 
Whereas the ratios for BR{t — 3/i) enter in to the allowed region by the experimental 
bound for large enough Mq and/or M1/2, the ratios for B{r —>■ fij) are well above their 
bound for all Mq and M1/2 values explored. The main point again is the particular value of 
62 with large |^2| and large arg(^^2); which generates large rates. 

With the purpose of exploring other choices of the mSUGRA parameters, we have also 
generated results for the specific value Aq = —100 and found very close predictions to the 
Aq = case, the lines in the plots being nearly undistinguisable respect to this case. We 
have also run the alternative case of sign(/i) < 0, and found again very close predictions to 
the sign(/i) > case, with the lines in the plots being undistinguisable from this case. 
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FIG. 14: Dependence of BR{t jj, ^ fi'^) and BR{t — > jj.'y) with Mo and M1/2 in scenario B 
with hierarchical heavy neutrinos and complex R, for [mN-i^^mN^iT^Nz) = (10^,2 x 10^, IQ-*^^) GeV, 
tan/3 = 50, |6'2| = 2.8 and arg(6'2) = 7r/4 (6*1 = $3 = 0). (a) Upper-left panel, BR{t 11-11-11+) 
as a function of Mo for M1/2 = 100 GeV, (b) upper-right panel, BR{t ^7) as a function of Mo 
for Ml/2 = 100 GeV, (c) middle-left panel, BR{t 11-^-11+) as a function of M1/2 for Mo = 200 
GeV, (d) middle-right panel, BR{t — /i7) as a function of M1/2 for Mq = 200 GeV, (e) lower-left 
panel, BR{t — fi-fi-fi+) as a function of Mq = M1/2, (f) lower-right panel, BR{t — fij) as a 
function of Mq = Mi/2- In all the plots we take ^0 = and sign(/i) > 0. The horizontal lines are 
the experimental bounds. 
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FIG. 15: Dependence of LFV r ^ 3/i with |02| in scenario B with hierarchical heavy neutrinos, 
for different niNi choices. Solid line is for {m]\f-^,m]\f^,mN-i) ~ (10^)2 x 10*^,10^^) GeV, dashed 
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GeV, ^0 = 0, sign(;u) > and arg(02) = vr/4. The horizontal line is the experimental bound. 
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FIG. 16: Dependence of LFV /i decays with ^13 in degrees in scenario B with hierarchical heavy 
neutrinos and R = 1, for {mN^,miyi2,mN.^) = (10*^,2 x 10^,10"'^'^) GeV, tan/3 = 50, j4o = and 
sign(^) > 0. The upper lines are for Mq = 250 GeV, = 150 GeV and the lower lines are for 

Mq = 400 GeV, Ml/2 = 300 GeV. The horizontal lines are the experimental bounds. 

Finally, we have also tried another input values for the heavy neutrino masses. The 
results for BR{t 3/i) are shown in fig. ^1 Here we compare the predictions for the three 
following set of values, (m^, , m^^ , m^Va ) = (10^2 x 10^ 10^^) GeV, (10^°, 2 x 10^°, 10^^) 
GeV and (10^,2 x 10^, 10^^) GeV. We conclude, that the relevant mass is the heaviest one, 
mATg, and the scaling with this mass is approximately as the scaling with the common mass 
m^v in the degenerate case. Because of this, the rates for the two first sets are nearly 
undistinguisable, and the rates for the third set are about four orders of magnitude below. 
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Last but not least, we consider the very interesting case where the angle ^13 of the Umns 
is non vanishing. It is known that the present neutrino data still allows for small values of 
this angle, 6*13 < 10°. The dependence of BR{^~ e~e~e~^) and BR{^ — > 67) with this 613 
is shown in fig. where we explore values in the < ^13 < 10° range. We choose these two 
channels because they are the most sensitive ones to this angle. For this study we assume 
the most conservative choice of i? = 1, and set the other parameters to the following values: 
tan/3 = 50, Aq = 0, sign(/x) > 0, and {mN„mN2,mNs) = (10^2 x 10^ 10^^) GeV. The 
upper lines are for Mq = 250 GeV, M1/2 = 150 GeV and the lower ones for Mq = 400 GeV, 
Ml/2 = 300 GeV. We conclude that, for this choice of parameters, values of ^13 larger than 
1 degree are totally excluded by the data on LFV /i decays. It is a quite stricking result. 

In summary, we obtain in the hierachical case much larger rates than in the degenerate 
one, and one must pay attention to these values, because the rates in several channels do 
get in conflict with the experimental bounds. More specifically, the choice of a complex 
R matrix with large modules and/or large arguments of 61 and/or 62 and a light SUSY 
spectrum is very constrained by data. We also confirm that the experimental upper bounds 
of the processes Ij li'~f are more restrictive than the Ij ^j^^i^t ones but all together 
will allow to extract large excluded regions of the mSUGRA and seesaw parameter space. 
A more precise conclusion on the excluded regions of this parameter space deserves a more 
devoted study. 

VI. CONCLUSIONS 

We have shown in this paper that the LFV r and fi decays do provide a very efficient 
tool to look for indirect SUSY signals. Whereas the predicted rates for these processes are 
negligible within the SM, the SUSY scenario considered here provides in contrast significant 
rates which are at the present experimental reach for some of the studied channels. This 
scenario consists of the well known mSUGRA extended with three right handed neutrinos 
and their SUSY partners, and with the needed neutrino masses being generated via the 
seesaw mechanism. The reason for these significant rates is because of the important lepton 
fiavor mixing that is generated in the slepton sector due to large Yukawa neutrino couplings, 
which is transmited via the RGE running from the large energies down the electroweak scale. 

With the motivation in mind of testing SUSY we have studied exhaustively in this work 
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the particular decays r — > 3//, r — > 3e and // — > 3e, and the correlated radiative decays 
T — > //7, T — > 67 and /i 67. All of these channels have quite challenging experimental 
bounds and they are expected to improve in the future . We have explored the dependence of 
the branching ratios for these LFV processes with the various parameters involved, namely, 
the mSUGRA and seesaw parameters. We have computed and analyzed in full detail all the 
contributions from the SUSY loops to the Ij ^i^i^t decays. Our analytical results for 
these decays correct and complete previous results in the literature. In particular we have 
presented the results for the separate contributions from the 7-penguin, the Z-penguin, the 
Higgs-penguin and the box diagrams and shown explicitely the 7-penguin dominance. In the 
numerical estimates we have presented results for both the I J — >■ 17^7 and the correlated 

Jill 

radiative decays Ij —>■ lij. 

For the degenerate heavy neutrinos case, we have got rates for all the studied LFV r and 
II decays that are below the present experimental upper bounds. The largest rates we get, 
within the explored range of the seesaw and mSUGRA parameter space, are for the r decays. 
Specifically, BR{t — > /j,^) ~ 10~^ and BR(t~ — > //"/x"//"*") ~ 3x 10~^^, corresponding to the 
extreme values of tan/3 = 50 and uin = 10^^ GeV and for the lowest values of Mq and M1/2 
explored. The case of hierarchical heavy neutrinos turns out to be much more interesting. 
First of all, we get much larger branching ratios than in the previous case and secondly they 
are in many cases above the present experimental bounds. 

We have analyzed in detail the behaviour of the branching ratios with the mSUGRA 
and seesaw parameters also in the hierarchical case. The largest ratios found are again for 
r ^ /^7 and t~ n^fi^^^ decays. All the LFV r and decay rates are mainly sensitive 
to tan/9, the heaviest neutrino mass rriN^, which we have set to 171^,^ = 10^^ GeV, and the 
complex angles in the R matrix 61 and 62, which have been taken in the range 3 < tan f3 < 50, 
< l^il < 3 and < arg(^j) < 7r/4. For the values of these parameters at the upper hmit of 
this studied interval we have found that some of the predicted branching ratios are clearly 
above the corresponding experimental upper bounds. The most restrictive channels being 
/X —>■ 67, fi ^ 3e and r fij. Therefore, we get in this region important restrictions on 
the posible values of the mSUGRA and seesaw parameters. In particular, for 62 = 2.86*'^/^, 
we get that the whole studied range of 100 GeV < Mo,Mi/2 < 800 GcV with tan/? = 50 
is totally excluded by r ^ /X7. Values of Mq and M1/2 in the low region below 250 GeV 
are also excluded by t" — > data. The case of 9i is even more restrictive, because 
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the predictions for fi — 67, — 3e and r — > /i7 totally exclude a light SUSY scenario, for 
practicaly all 61 values. 

Perhaps, the most striking result is that even for the most conservative choice of i? = 1, 
that is 61 = 62 = 63 = 0, there are also important restrictions at low Mq, M1/2 and large 
tan f3 values. In particular, for tan (3 = 50, values lower or equal than Mq = 250 GeV and 
Ml/2 = 150 GeV are totally excluded by r — > ^•y, — 67 and /i~ e~e~e~^ data. 

For this conservative choice of i? = 1 we have also found the surprising result that both 
/i — s> 67 and /i~ — > e^e~e~^ place important restrictions on the allowed values for the Umns 
angle ^13. For values lower or equal than Mq = 250 GeV and M1/2 = 150 GeV and for 
tan/5 = 50 and m^vg = 10^^ GeV, we get that values of 6*13 larger than 1 degree are not 
allowed by these LFV data. 

In conclusion, it is clear from these results that the LFV r and fi decays studied here do 
restrict significantly the mSUGRA and seesaw parameter space. A more refined analysis of 
the restrictions on this multidimensional parameter space, deserves a further study. 
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APPENDIX A 

In this appendix we collect the Feynman rules for the interactions that are relevant in 
this work. They are expressed in the physical eigenstate basis, for all the MSSM sectors 
involved: sleptons Ix {X = 1, .., 6), sneutrinos ux {X = 1, 2, 3), neutralinos Xa = 1? 
charginos Xa ~ 1' 2) and the neutral Higgs bosons Hp {p = 1, 2, 3) = H^, A^. 
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1. Photon interactions 



The Feynman rules for the photon interactions that are used in this work are given by, 




'ie {Pf, + g^) SxY 



2. Neutralino interactions 

The Feynman rules for neutralinos that take part in the one-loop diagrams computed 
here are the following: 
where 

N!ax = -9V2 I ^M^cos/j ^^^^S^-^)^ + tan^H/iVl.i^g, ,)^| (Al) 

(A2) 



C is the charge conjugation matrix and Pl,r = ^^2^- Here i?*^') and are the rotation 
matrices in the charg 
be found in 3- 



matrices in the charge slepton and neutralino sectors, respectively. The definition of A^ can 

nn 
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i {NIaxPl + NI'axPr) 




fx 




^ {N!axPl + N^axPr) C 



i-x 




x\ 



3. Chargino interactions 



The Feynman rules for the chargino interactions are given by 



yx 




Xa 



{^iaxPl + C^axPr) 



T^X 




Xa 



V-'iAX^L + ^iAX^R) 



where 



AX ~ y ^ U A'7-rL 



V2Mw COS p 



^2-'i'(l,2,3)X 



C[ax - -9yAiR\i'2, 



(1,2,3)X 



(A3) 
(A4) 



and R'-''\ U and V are the rotation matrices in the sneutrino and chargino sectors, respec 



tively. The definitions of U and V can be found in j21l . 



221. 
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4. Z boson interactions 



The Feynman rules for Z boson interactions are given by, 



where 



where 



where 




E 



L{n) 
AB 



9 



E 



R{n) 
AB 



COS 9w 
9 




liL 




AB — 


- 




cw \ 


IIR 


9_ 


AB — 






Cw 



^A3^^B3 



E 



L(c) 
AB 



E 



m 

AB 



9 



COS 9w 
9 



^AB 



COS 9 



'^AB 



w 



_9_ 

Cw 

_9_ 

Cw 



2~ ^w] Ua2Ub2 - c^^UXiUbi 



Ixip) ^ 
Iriq) 



KVn + <ln)QxY 



--E 



, ^2 \ d(0* d(0 , ^2 p(0* p(0 

- "T •'VK / ^2k-l,X^2k-l,Y ''W^2k,X^2k,Y 
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XY 




where 



where 



9 



2c 



)XY 



W 



7{l) 



Cw 

Cw 



(AlO) 

(All) 
(A12) 



We have used here and everywhere the short notation sw — sin^^K and cw — cos9w 



5. Higgs boson interactions 



The Feynman rules for the three neutral Higgs bosons read as, 




where 



and 



D 
D 



ip) 

L,AB 
ip) 

R,AB 



9 



sin P 
9 

sin (3 



'^BA^! 



(P) 
5 



+ ^, , (^2 ^BA 



2Mi 



w 



^ba(^5 - ^BA'^2 + 2Mw 



Qab = 2 ~ tanOwNBi) + iV^s (iV^s - tan^^yiV^i)] 



^AB — 



2M, 



[M;Na2Nb2 + M*NaiNbi - II* {NasNb4 + Na^Nbs)] 



(A13) 
(A14) 

(A15) 
(A16) 
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where 



and 
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sin /3 






sin/3 



m. 



a 2 Ob A 



QabC^^* — RaB(7^2^* + 



2Mw 
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2Mi 
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V2 
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Ua2Vi 



Bl 



2M, 
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[M;UAlVBl+fJ,*UA2VB2] 



lx{p) 



,p(i) 

XY 



where 



^XY — 



with 



" (P) R*(0 r{0 _4_ Ap) D*(0 d(0 , ^(P) R*{0 d{0 , „(P) D*(0 D 

yLL,e^lX ^lY ~r 9 B.R,e'^2X ^2Y yLR.e^Hx ^2Y y_RL,e-"'2X 

, „(P) o{0 , ^(P) n{0 , JP) j^*H) n{0 , ^(p) p*(0 p(0 

"T yLL,n'^^3X ^^3Y "T yRR,fi'^Hx ^HY T" yLR,n^Hx ^HY yRL,^i'^AX ^^Y 

, ^(p) p*(0 p(0 , ^(p) p*(0 p(0 , ^(p) p*(0 p(0 , ^(p) p*(0 p(0 

"I" yLL,T^bX -"-SY "I" yRR,T^6X ^GY ' yLR,T^5X ^6Y ' yRL,T^eX ^5Y 



(P) 
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(p) 
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cos ^ 



2 ^(p) 



- - sm^ dw + 77 — — 7 
2 y Mw. cos 



^(P) 



cos fvi^ 



I -(7l A; - (72 /XI 



Mw cos 



^(p) 
(74 



2Miy cos /3 



(p)* 

9lr,i 



with A/ = (^;)"/ {Yiy^, i = 1, 2, 3 for Z = e, /x, r, respectively. 
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iG' 



XY 



where 



with 



G 



XY — 9 



(p) (p) dC*^) _i_ d*W dC^^) 

yLL,v-^lX ^lY "T yLL,v^X ^2Y "r yLL,v^SX ^SY 



ip) 



2 cos 



(A26) 
(A27) 




where 



q(p) 



(p)* 



(j>)* 



(A28) 
(A29) 
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In all the above equations, 



» 



Go = 



(J4 



^ sin a ^ 



— cos a 



^ ism. (5 y 
^ cos a ^ 



sm a 
—i cos /? 

sin (a + /9) 

— cos (a + /?) 



— sin a ^ 

cosa 



V 







— cos {P — a) 
sin (/3 — a) 
i cos 2/3 



/ 



(A30) 



(A31) 



(A32) 



(ASS) 



(AS4) 



and Hp{p — 1,2,S) = We have also used here the standard notation for the 

MSSM soft-gaugino-mass parameters Mi^2 and the parameter. 



APPENDIX B 

In this appendix we present the analytical expressions of the loop-functions for the cal- 
culations of the I J I'i^l^l'j^ decays. In these expressions we neglect the external fermion 
momenta/masses which for the present computation works extremely well. That is, 

B{k'^,ml,ml) ~ B{0,ml,ml) = B{ml,ml) (Bl) 

C{kl,kl,ml,ml,ml) ~ C {0,0, ml, ml, ml) ^ C {ml, ml, ml) (B2) 
D{kl,kl,kl,ml,ml,ml,ml) ~ D{0,0,0,ml,ml,ml,ml) 

= D{ml,ml,ml,ml) (BS) 
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1. Two-points functions 



The analj^ical expressions for Bo and Bi functions are the following: 



,2 ^2 I ^2 



Bo{ml,ml) = -logm^ + ^ (^4) 

777,2 — mf 

1 1 m? — mo + 2m?log ( ^ ) 

i^iK,rn^) = + llogml -— (B5) 

^ ^ 4 (mf — 777.2) 



2. Three-points functions 



The expressions for the three-points functions used in this work are given by, 

1 / 777^ log 777^ — 1712 ^2 ^1 ^og ni'l — ml log ml 



/ 2 2 2\ 

Co{mi,m2,m^) 

ri I 2 2 2\ 
60(7771,7772,7773) 

611(7771,7772,7773) 



X 



X 



2 2 
7772 ~ "^3 
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2 2 
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2 2 
777i — 7773 



7772 ~ "^3 



777i log — nil log /;i2 777i log llll — 111^^ log m\ 



2 2 
777J — 7772 



2 2 
777i — 7773 



(B6) 
(B7) 



m. 



2{ml — ml)'^{ml — mf)'^{ml — mf) 



— (ml — ml) (ml — 7773) (7773 — m,l) + mlm^(2w^ — m^) log — | 



777n 
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777n 
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3 J 



(B8) 



612(7771,7772,7773) = 



2(777i — ml)(ml — mlY(ml — m^)^ 



mf, 



(mi — 7772) (7r7i — 7773) (7772 — m^)m^ + 7^737^73(27771 — 7773) log — ^ 



mi 



A I 4 , 777 1 2/ 2 2\i 1 1 t-i 

+ 777i I 7772 ^^g — 2 + m^[—2m2 + 7773) log — ^ 



ml 



777? 



7775 



624(7771,7772,7773) = -60(0,0,7771,7772,7773) 



(B9) 
(BIO) 
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3. Four-points functions 



Finally, the four-points functions have the following expressions, 

mf log ml m^ log m'^ 



Do[m^,m2, m^, m^ 



f\ / 2 2 2 2 

JjQ[mi,m2, m^, m^ 
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(B12) 
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